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Abstract 



Q . We propose a new formulation of the probability distribution function of wind driven water 

I surface slope with an a-stable distribution probability. The mathematical formulation of the 
(N . 

' probability distribution function is given under an integral formulation. Application to represent 

Q_i' the probability of time slope data from laboratory experiments is carried out with satisfactory 

O , results. We compare also the a-stable model of the water surface slopes with the Gram-Charlier 

• ■ n 

^ I development and the non-Gaussian model of Liu et al[15!]. Discussions and conclusions are 

• 1— I ' 

c/3 ' conducted on the basis of the data fit results and the model analysis comparison. 

^ : 

I Keywords: Wind driven water surface wave slopes. Probability density function. Non- 

J> ' Gaussian distribution. Alpha-stable distribution. 

^ : 

■ During the last three decades, remote sensing of ocean surface has been extensively 

! investigated to study physical processes on the sea surface and air sea interaction 

o ■ 

I ' mechanism. One of the most important quantities on which depends the ocean remote 

00 ; 

O ■ sensing process is the surface wave slope probability distribution that is known to be 

> '■ 

• ^ . non-Gausssian. In the past the general way to describe such probability distribution 

X: 

H was to use an expansion technique of the Gaussian law known as Gram-Charlier 

series. In the present paper, we propose a new formulation of the surface wave slope 
using alpha stable probability distribution. Application to data from laboratory 
experiments and comparison of the alpha stable slope model with existing models 
are carried out. Our results show the potentiality of the alpha stable laws to describe 
the non-Gaussian variability of water surface wave slope. 



Electronic address: maminirina.joelson@univ-avignon.fr 



1 



I. INTRODUCTION: 



A knowledge of the sea surface slope distribution and its related statistical quantities is es- 
sential to study the sea surface using remote sensing processes. The use of radar backscatter 
to characterize the sea surface has been widely investigated experimentally and theoretically. 
At nearly normal radar beam incidence, the importance of the slope distribution is clearly 
exhibited by the simplest model of cross section radar backscatter based on geometrical 
optics approach formulated by Valenzuela[27j] as: 




where 6 is the radar incidence angle, |i?(0)p is the Fresnel reflection coefficient for normal 
incidence and f{Cx,Cy) is the probability distribution function (PDF) of the sea surface 
slope at specular point {(x,Cy)- Determination of the surface slope distribution has been 
carried out in a large part with semi-empirical models derived from data observations. The 
commonly used approach for such purpose was the truncated Gram-Charlier series under 
the assumptions of the near-Gaussian case and the weakly nonlinear interaction proposed 
by Longuet-Higgins [16 ] . The Gram-Charlier series was attempted to include two additional 
factors: the skewness and the peakedness that are respectively due to second and third 
order wave-wave interactions. Despite the multiple adaptations made on the Gram-Charlier 
model, it is now well-known that this approach fails to work in the range of large slope or 
large surface elevation. In other word, this model is not suitable to represent the reality 
when large events have higher occurrence. To improve the results from Gram-Charlier series, 
several works were made by different authors. Among the existing models, those of Liu et 
alfisl. differs by its formulation and its mathematical nature. Indeed, Liu model does not 
appeal to an expansion technique and is well suitable to represent extreme events. However, 
this model is originally a symmetrical law including the Gaussian probability and the Cauchy 
law as limit cases. Owing to the symmetry property, the Liu model is not able to reproduce 
the skewness effect. Liu extended empirically the model for this purpose. 

In the present paper we propose to use the a-stable law to describe the PDF of surface 
slope. The stable law depends on four parameters: the index of stability, the skewness, the 
scale and the location parameters that allow one to represent a probability law of random 
variations including asymmetry, peakedness and peak shift. Also, the a-stable law includes 
as limit cases the Gaussian probability and the Cauchy law. Generally speaking, stable law 




aoi9) = 7rsec'ef{CxXy)\Rm' 



(1) 
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belongs to the family of heavy tails probability laws that is especially well adapted to take in 
account the large events of random variations. Such strong random variations can be found 
in fields or laboratory observations of the water surface roughness (Onorato et alj2o|). 

Evidences from in situ and laboratory experiments of non-Gaussian heavy tails distribu- 
tion of wind driven water surface waves and their slope may be found in many articles in 
the literature. The first argument on the non-Gaussian behavior of wind driven water sur- 
face waves and their slopes concerns the existence of asymmetry prop erty of wind wavefield 
that can be quantified by higher order spectral analysis (Leykin et al 13|,Joelson et al 10|). 
Another hallmark of this non-Gaussian behavior is based on the power law observed in the 
spectrum of wind driven water surface waves. We cite here some works emphasizing this 

n 

fact. Experiment observations was conducted by Huang[8| on time series of sea surface that 
suggest a random dynamic as a spectral power law associated to fractal description. Mellen 



and al 



171] performed the study of random nonlinear dynamics of laboratory wind waves by 



the Levy index and the co-dimension parameter that quantify the multifractal property of 
the water surface. As pointed out by Solow {2^, power law behaviour can arise from random 
distributions having heavy upper tail. A fundamental consequence is that random rough 
surfaces with slowly decaying power spectral density can have infinite slope variance (War- 
nick et al[28]). Such case of slope behavior is clearly beyond the Gaussianity assumption and 
is only reached by the family of heavy tail probability distribution model. The a-stable law 
belongs to this class of model. Heavy tail distributions of wind driven water surface slope 
obtained from experiment observations can be found in different works (see for example Cox 
and Munk j^, Lui et alfisj], Chapron et alQ]). 

Direct simulation methods such as a Monte-Carlo simulation based on nonlinear hydro- 
dynamical model are able to reproduce pseudo-random samples that present similar aspects. 
Recently, Soriano et al[25| studied the effect of this non-Gaussian character on the radar 
Doppler spectrum of the random surface at grazing angle and concluded on the importance 
of the non-Gaussian character. Despite these results and the fact that stable laws have been 
discovered several decades ago, there are only few applications concerning the a-stable dis- 
tribution of sea surface. In particular, we cite here the work of Guerin[3] on the scattering 
on random rough surface with a-stable law in which he gave an extension of the expression 
([1]) for the a-stable surface roughness. 

The present study, in our opinion, is the first attempt to apply the stable probability 
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law on random water surface slope representation. The paper is organized as follow: in §2, 
we recall the main properties of the a-stable law. In §3, the application of stable law to 
fit instantaneous slope of wind driven water surface is presented. The §4 is devoted to the 
comparison of the stable PDF with the Gram-Charlier model and the Liu model. In §5, we 
present our conclusions on the results of the study. 



II. STABLE PROBABILITY DISTRIBUTION 



The stable laws were introduced by Paul Levy ij] during his investigations of the behav- 
ior of the sums of independent identically distributed random variables. Then the theory 
on stable laws was developed by many authors and laid out in many books and articles. We 
recall here some basic results without proof. A stable distribution law is determined by four 
parameters: an index of stability a, a location parameter fi, a skewness parameter j3 and a 
scale parameter 7. The most common formulation of such law is given by the characteristic 
function which represents the inverse Fourier transform of the probability distribution func- 
tion. A random variable X is said to have a stable distribution (noted as X ~ S'a(/3, 7, //)) 
if its characteristic function satisfies 



^x{t) = exp {ifit - 7"|t|" [1 - i(3 sign(t) W(a, t)]) 



(2) 



where 



and 



W(a,t) 



tan(^) if a ^ 1 
--log \t\ if a = 1 



1 if t > 
sign(t) = I if t = 
-1 if t < 

The stability index also called the characteristic exponent or the tail index satisfies < a < 2 
and the skewness parameter — 1 < /3 < 1. The location parameter fi is defined only when 
the stability index is a > 1 (/i G R). In this case, the first absolute moment is given by 
< \X\ >= ^r(l - 1/a) where r(.) is the gamma function. If /5 = 0, then the distribution 
is symmetric around fi. The scale parameter also known as dispersion parameter has a 
real positive value 7 > that determines the width of the distribution law. The location 
parameter /x describes the shift of the peak. The tail index a determines the rate at which 



the tails of the distribution taper off. When a < 2, the variance of the random variable X 
is infinite and the tails exhibit a power-law behavior. More precisely, using the central limit 
theorem, it can be shown that: 

' lim rPr(X > t) = CJl + 

< (3) 

lim rPr(X < ~t) = Cail - /?)7" 

^ t — *oo 

where = ^r(Q;)sin(^). It follows from Q that in general, the p^^ order moment of a 
stable random variable is finite if and if only p < a. 



Probability density function 

There are only three cases of stable distributions for which a closed-form expression of the 
PDF is known. This lack of closed form formulas for the PDF has a negative consequence on 
the popularity of the stable distribution. By applying Fourier inversion to ([2]), an integral 
form of the PDF is obtained as: 



fa,l3, 



7.M 



X 



exp(— 7"z") cos 



exp(— 72;) cos 



(x — fi)z + Pz" tan(- 



dz 



if a ^ 1 



IT 



[x- ij,)z + f3z- \og{\z\) 
if a 



dz 
= 1 



(4) 



Expression (jlj) is known to have analytical solution for the following cases: 
-for a = 2 corresponding to the Gaussian distribution. 



f2{x) 



1 



27r7 



exp 



272 



(5) 



in this case, the skewness parameter is irrelevant. The variable X ~ 5*2(0, 7, /x) reduces to 
the Gaussian random variable X ~ A/'(/i, V^'j) 

-for a = 1 and P = which leads to the Cauchy distribution 
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TT 7^ + (a; — /i)2 

The Cauchy random variable writes as X ~ 6*1(0, 7, /x) 
-for a = 0.5 and [3 = 1 giving the Levy distribution 

exp{-i:} 



(6) 



/l/2(x) = y^ 



X 



3/2 



(7) 



For the general cases, numerical approximations of (jl]) was widely used by several authors 
to represent stable PDF. Different methods exist for this purpose with various level of 
accuracy^ In this paper, we adopt the direct integration method developed by Nolan|l9| and 



Weron 



301 1 using an other formulation of the expression (jl]) called the Zolotarev's formulas 
3l| . The main interest of the Zolotarev's formulas is that instead of expression (jlj), they 



do not include infinite integral and then are well adapted to numerical computations. In 
which follow, we give the Zolotarev's formulas of the PDF of a stable random variable 
X ~ S'q(/5, 7, yu). First, we consider the case of a random variable Xq ~ S'a(/5, 1,0). By 
setting ( = —(3 tan^, the PDF of Xq can be expressed as: 
if a 7^ 1 



V{e,a,(3) exp{-(x-C)^ V{e,a, (3)} de, 

if X > C 
if X = ( 

71 

fa -f3, l,oi-x), if X < C 



7r\a — 1\ 

r(i + cos(0 



(8) 



if a = 1 



/iAi,o(a;) 



!|/3| ^-^P^ 2/3^ 



I 7r(l + x2) 



Vie, 1,(3) exp(- exp(-— ) V(^, 1, 

if /3 7^0 
if (3 = 



where ^ = ^ arctan(— (^) and 
V{9,a,(3) = 



cos ( ■ T/, 

Vsm a[t, + 



-1 cos {a^ + (a - 1)^) 



cos 



2 /f + 



TT 



a ^ 1 
1 



a 



From the probability distribution given by ([8]), one can easily build the PDF of the 
random variable X with the help of the stability property, X = 7X0 + yU ~ Sa{P, 7, /x). 

Parameters estimation 
The first step on modelling data with a stable law requires the estimation of the four pa- 
rameters that determine the distribution properties. The estimation process suffers from the 
lack of known closed form of the PDF. As a consequence, standard methods based on maxi- 
mum likelihood principle would not be efficient with regards to the stable PDF dependance. 



However, there are other numerical ways that do not depend directly on the PDF knowledge 
to evaluate stable parameters like the sample quantile methods or the sample characteristic 
function methods. A comparative analysis and discussions about these methods may be 
found in Borak and al[l]. Among these existing methods of estimation, we will use the 
sample characteristic function developed by Koutrovelis 12(] and performed by Kogon and 
Williams The method is based on regression type on the log-characteristic function. In- 
deed, the logarithm of the real part and the imaginary part of the log-characteristic function 
are linear and then, give rise to a regression model with the data. 

log (-Re[^(t)]) = a \og\t\ + a log7 

(9) 

Im[*(t)] = -^t- [3 sign(t) W(a,t) 

where \l/ is the log-characteristic function of X and W(.,.) is defined as in expression (21). A 
more complete description of the method and an extended study may be found in [llj. 

In the following section, the stable PDF model will be applied on slope data from labo- 
ratory experiments on wind driven surface water waves. 



III. MEASUREMENTS AND DATA ANALYSIS 



Experiments 

The data of interest are issued from experiments made in the wind wave facility of the 
IRPHE-IOA Laboratory. Basically, the facility consists of a water tank of 40 m length, 3 m 
width and 1 m depth. The facility can be described as a combination of a wind tunnel with 
a wind-wave tank. A schematic representation of the air-sea interaction simulation facility 
is given in FigUl The facility is equipped with a submerged wavemaker that allows to 
produce mechanical waves from capillary ripples to breaking waves. However in the present 
work, we focus only on wind wave fields. By means of an axial fan, air flow with velocities 
ranging from m/s to about 15 m/s generates wind wave fields. Different types of surface 
waves can be produced by the action of the air flow. Noteworthy that the facility includes 
particular devices intended to control various parasitical effects which could coexist with 
the physical process under consideration. In particular, a permeable wave absorber is set 
up at the upwind end of the water tank to prevent the wave reflection. Within the range 
of wave frequency of interest in the experiments on surface waves, the reflection coefficient 

7 



is estimated to be insignificant. As far as wind generated waves are concerned, a specific 
design lias been adopted for a smootli joining of tlie air flow and the water surface. The 
design prevents air flow-separation at the facility entrance test section. This insures a natural 
development of the turbulent boundary layer over the water surface. In addition, later quays 
disposed in the water tank favor the three-dimensional evolution of the wind waves. Studies 
with more details about these specific devices are given in a number of publications (see e.g. 
Coantic et al [3]). We concentrate here on the time evolution of the water surface deflection 
level at a given position. Measurements of water surface height ri{t) were performed by two 
capacitance wave gauges separated by a 5cm along wind direction. These probes using a 
capacitance wire gauge of 0.3 mm outer diameter were set up on a carriage moving along the 
wave tank. Experiments were carried out under various conditions, whose main parameters 
are the wind velocity and the distance between the gauge and the facility entrance section 
(hereafter denoted fetch). Such distance determines the length action of the air flow over 
the water surface and is well known to be of crucial importance on the rate growth of surface 
waves and on its nonlinearity. The time derivative d?7/dt of the water surface heights were 
obtained by using analog derivators. 




FIG. 1: A schematic representation of the large air-sea interactions facility. 

The corresponding processed series are made of n = 36000 samples gathered at a con- 
tinuous rate of 200 Hz. The water surface wavefields of interest here are known to have 
temporal as well as spatial dynamics in which dispersion, nonlinearity, dissipation (by kine- 
matic viscosity or by breaking) and wind action are involved. The results presented herein 
are limited to some aspects of the temporal evolution at the given space locations. Thus, 
instantaneous slope s(t) on the upwind direction of the water surface waves could be de- 
rived as s(t)=-l/c drj/dt where c corresponds to the phase velocity of the dominant wave 
frequency. This slope measuring method was tested by comparison with direct measurement 
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Case Fetch (m) U lo (m/s) a (3 j fi 



I 



2 



6 



1.7219 -0.0816 0.5434 0.0056 



II- 1 6 



6 



1.8688 -1.0000 0.6384 -0.0869 



II-2 13 



3 



1.9822 -1.0000 0.6972 -0.0117 



III 26 



13 



1.9178 -0.9488 0.6481 -0.0051 



TABLE I: Estimation of stable parameters (a, /?, 7, ^u) from the data series. 



performed by a laser slope gauge (see Jhane and Riemer A good agreement was found 
particularly for capillary and gravity waves at moderate wind speed. 
Data analysis 

In this section, we estimate the PDF of time water surface slopes at different wind speed 
for fetch varying from 2m to 26m. Note that wind speed is estimated at 10m height above 
the surface(Uio). Among the large possibilities of combinations from wind speed and fetch 
values, we present here some cases that are indicative of the typical behavior of water surface 
slope. In which follows, we focus on slope behavior of three cases: capillary waves, capillary 
gravity waves at moderate fetch and low wind speed and gravity waves at large fetch and 
high wind speed that are typical of fully developed sea surface. In order to compare the 
different runs of study, the data series are normalized with theirs rms and mean values. The 
results of the fit are summarized in TabO 

The first line of Tab H] corresponds to the case I of capillary wind wave slopes at moderate 
wind speed and at a very short fetch. In this case, we find a stability index of 1.72 that is 
indicative of a strong non-Gaussian character. A very weak value of the skewness parameter 
is also found traducing a lack of asymmetry on the waveform. The scale parameter is at a rel- 
ative low value compared to the normalized Gaussian standard deviation(~ «S'2(0, 0)). 
This may be related to a nonlinear effect involved on the generation of capillary waves by 
wind. The characteristics of the capillary wave slopes are shown in Figj2l The FiglJla de- 
picts a sample of the capillary wave profile and its slope while the corresponding probability 
distribution is shown in the part (b) of the figure with a normalized Gaussian law as a ref- 
erence. At larger fetch value corresponding to case II, the stability index grows. The scale 
parameter also increases. A striking feature is that the skewness parameter reaches to its 
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FIG. 2: Fetch=2m. Uio = 6m/s. Part(a): Thick line: Normalized sample of surface wave. Thin 
line: Normalized instantaneous slope sample. Part(b): 'o': PDF of slope experiment data. Solid 
line: a stabe PDF. Dashed-dotted line: Normalized Gaussian PDF 

negative limit value. This is in good agreement with time evolution of the slope presented 
in FigS.a. Observing also that the corresponding surface waves are dominated by a Stokes 
component, we suggest the possible role of the Stokes effect in the asymmetric behavior of 
the slope and in the limit value obtained for the skewness parameter in the present case. 
However, in the general case, relation between asymmetry in the surface and in the slope 
variation is not easy to establish especially for a random wavefield. The Stokes effect yields a 
waveform with steeper crests increasing locally slopes and flatter troughs. Results of run at 
fetch of 6m and wind speed of 6m/s (case II-l) are shown in FigJSl The asymmetry property 
is clearly visible on the slope samples (Fig|3la) and is well represented by the PDF model 
(FiglHlb). Modulational instabilities begin to occur in the wavefield owing to the nonlinear 
wave- wave interaction processes as the Benjamin- Feir instabilities. The sample of the wave 
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signal presents wave trains formation that are reminiscent of these instabihties effects (see 
FigEla). Such nonhnear effects are known in this case to dominate the effects of randomness. 
This would explain the reduction of slope magnitude and also the increase of the stability 
index to the near-Gaussian situation. By increasing the fetch value but taking a low wind 
speed (case II-2), the rate of input wind energy will be weaker than the rate of the self- 
nonlinearity of the wavefield. Thus, one can better put forward the wave modulation effect 
on the slope distribution and on its PDF. The corresponding PDF is shown in FigJUb with a 
very similar form to the normalized Gaussian except for the negative tail that is indicative of 
the presence of Stokes effect. Finally, at high wind speed and large fetch corresponding to 
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FIG. 3: Fetch=6m. Uio = 6m/s. Part(a): Thick line: Normalized sample of surface wave. Thin 
line: Normalized instantaneous slope sample. Part(b): 'o': PDF of slope experiment data. Solid 
line: a stabe PDF. Dashed-dotted line: Normalized Gaussian PDF 

the case III, one can surprisingly find that values of the index stability but also the skewness 
and the scale parameter reduce reaching again a clearly non-Gaussian situation. In agree- 
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FIG. 4: Fetch=13m. Uio = 3m/s. Part(a): Thick line: Normalized sample of surface wave. Thin 
line: Normalized instantaneous slope sample. Part(b): 'o': PDF of slope experiment data. Solid 
line: a stabe PDF. Dashed-dotted line: Normalized Gaussian PDF 

ment with these values, from the FigOa, it may be found that the magnitude of the wave 
slopes grows locally at the crests and the wave signal behaves strongly asymmetrical. From 



the physical background, it was shown by Trulsen and Dysthe[26| that such situation may 
be explained from a theoretical model based on a modified nonlinear Schrodinger equation. 
It allows one to predict the stabilization of the Stokes effects by the strong wind leading 
to the suppression of the modulational instabilities. The PDF of the case III are given in 
Figinib in which it may be found that the stable PDF model reproduce well the distribution 
of the negative slope values and also the skewness of the slopes. However, at large positive 
values, the model fails to represent the data distribution. No clear interpretation was found 
for this difference but we believe that additional physical processes such as breaking wave 
would be involved. 
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FIG. 5: Fetch=26m. Uio = 13m/s. Part(a): Thick line: Normalized sample of surface wave. Thin 
line: Normalized instantaneous slope sample. Part(b): 'o': PDF of slope experiment data. Solid 
line: a stabe PDF. Dashed-dotted line: Normalized Gaussian PDF 
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IV. COMPARISON WITH OTHER MODELS OF SLOPE DISTRIBUTION 



In this section, we compare the properties of the stable law with the Gram Charlier and 
the Liu models. The Gram-Charlier model based on an expansion technique called Edgworth 
series was widely used in the past and permitted to develop a large part of the understanding 
of the non-Gaussian behavior of the sea surface. However, the limit of this model was earlier 
found by several authors. Many efforts were conducted to develop alternative ways from 
the statistical expansion technique. The Liu model was built in this context to improve the 
statistical representation of the sea surface slopes. 

The Gram-Charlier series 

The traditional method of interpreting the wind wave slopes has been to fit them to 
a Gram-Charlier series, which can be derived by assuming that surface waves are weakly 
nonlinear (Longuet-Higgins[16]). Basically, in this procedure the PDF of the wave slopes 



s(t) writes as: 



where 



1 + — Hgfs) + — RAs] 
6 -^^^ 24 ^ ^ 



(10) 



G(s) 



exp 



< s > 



V2i 



a 



th 

and a is the rms of the slope, when H„ is the n Hermite polynomial. The coefficients Cr, 
can be calculated from the statistical moments of 0gc(s) and G(s) as: 



cr" 



{s)-G{s)]ds 



Clearly, the coefficients C3 and C4 in expression ffTOj) are intended to include the skewness and 
the kurtosis that are the consequences of nonlinearities in the underlying physical processes. 
Of interest is that the PDF model becomes non-Gaussian. Note also that expression ( |T0|) is 
not always guaranteed to be positive, and is therefore not a valid probability distribution. 
The Gram-Charlier series diverges in many cases of interest. It converges only if 0gc(s) 
falls off faster than exp(— at infinity (Cramer|^). For water waves applications, the 
skewness and kurtosis coefficients (C3 and C4) are found to be related to physical parameters 
in particular to the air flow speed value at 10m altitude (noted as the Uio). Here we adopt 
the values of these coefficients used by Cox & Munk ^ and Plant 21 1 from a range of wind 
speed varying from 0.72 to 10.2 m/s (Note that jJ] uses the wind speed at height of 12.5m 
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instead of the traditional Uio). We fit tlie Gram-Charlier model with stable law using the 
Nelder-Meadjisl simplex method. The results of the fit are given in TabllTl The method is 
based on the minimization of the euclidian norm of the error e = 1 1 ia,i3,-y,fM ~ 0gc 

1 1 for 1000 

points varying from -5 to +5. From each case, values of the best stable parameters and the 
error e are shown. 





CM CM 

(m/sj a U3 U4 


a (3 


7 ^J' 


e 




0.72 


0.0005 0.101 0.127 


1.951 -0.330 0.695 -0.05 


le-i 




3.93 


0.0098 0.003 0.129 


1.957 0.022 


0.695 


1.9e- 


'2 


4.92 


0.0174 -0.080 -0.019 1.980 0.900 


0.708 0.040 


8.6e- 


-2 


6.30 


0.0170 -0.143 0.101 


1.974 0.880 


0.697 0.070 1.56-^ 


-1 


8.00 


0.0191 -0.156 0.173 


1.930 0.260 


0.690 0.07 


1.56-- 


^1 


10.2 


0.0357 -0.283 0.128 


1.910 0.400 


0.690 0.13 


2.7e- 


-1 



TABLE II: Estimation of stable parameters (a, f3, 7, /i) from the Gram-Charlier series used by Cox 
and Munk (1954) model 



As a general remark for all cases, the values of the stability index a are close to 2 that 
traduces the near-Gaussian assumption adopted in the Gram-Charlier model. The error e 
increases with the wind speed values. Rather than a simple difference from the calculation, 
these error values reveal a more profound departure on the nature of both models. More 
precisely, at low wind speed, when the occurrence of large values of slope would be rare, a 
better agreement is found in particular around the peak of the PDF. At high wind speed, 
the results of the Gram-Charlier series differ clearly from those of the stable law both at the 
tails of the probability and at its middle part. The Gram-Charlier model underestimates 
the extreme events giving a less heavy tails on the probability law (see Fig 16]). At the 
middle part of probability law, it overestimates the effect of the skewness. Note that no 
clear tendency is found between the stable parameter (3 and the skewness C3 . Finally, the 
location parameter fi is globally around zero and the scale parameter 7 is found at a quasi 
constant value around They are directly related to the normalized Gaussian law of the 
function G(s) included in the expression (fTOj) . 

The Liu model of slope PDF 

A more advanced model is developed by Liu et al[l5|] to fit the joint probability of the 
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FIG. 6: Results of the fit for U12.5 =8m/s. The sohd hne is the Gram-Charher PDF. Dashed line 
represents the stable law of the fit. A fixed value of the probability is drawn with dotted line that 
allows one to visualize the asymmetry of both laws. 

surface wave slopes. Liu model is also based on the weakly nonlinear assumption of the 
surface [l^ . The basic form of Liu slope PDF writes as: 



n 



2n(n - 1] 



1 + 



X 



+ 



-(n+2)/2 



(11) 



{n-1) ' {n- 1) 

where x,y are normalized slope components and n an integer parameter describing the 
peakedness of the surface slopes. Note that (ITT!) refers to the case of isotropic sea surface 
and is not always suitable for anisotropic real surface. An empirical formulation inspired 



by |4l[l6| is used by 



15j to handle the skewness effect. According to 15j, the expression 



f ITT]) is an improvement of the Gram-Charlier distribution in particular for the range of large 
slopes. For convenience, here we consider the expression for y=0, i.e. we focus only on 
upwind direction. As a first remark, the positivity of the expression (fTTj) is obviously guar- 
anteed however it may be found that the expression of f2{x,0) does not represent a valid 
probability density function for any value of n. Indeed, the integral of this expression is not 
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equal to unity. We addres this point by adopting a normalized form of the Liu PDF as: 

,2 n(-f+i)r(^-i) r-T 



1 



1+ " 



2 ''V m-1 
■p/m 3\ 



[m- I] 

The expression (1121) is defined for m > 4. For m=4, one can easily find that 



(12) 



Expression (|T3l) belongs to the Cauchy stable law with X ~ 5'i(0, v^, 0). On the other hand, 
taking the logarithm of (fT2!) . the limit for large value of m writes as: 

log C(^) ^™^oo log [l + ^) (14) 

which reaches by simple series expansion to the normal distribution of X ~ A/'(0, 1) or also 
the stable law with X ~ 5*2(0, -^,0). The limit cases of the Liu formulation belong to the 
symmetrical stable distribution. 

For intermediate values, the Liu PDF differs slightly from the stable distribution. In order 
to illustrate this relation, we compute a numerical adjustment of the stable parameters a 
and 7 to the expression f|T2|) for m=4.5. The adjustment is based on the minimization of the 
euclidian norm || fa,o,7,o(2;) ~ 0l (^)ll — ^ 3000 points from x=-15 to 15. The results are 
shown in FigE It may be found that the major difference between the two laws are clearly 
localized at the medium range values of the probability structure. 

By varying the peakedness parameter m, this difference is retained systematically. An- 
alytically, one can also decompose the modified Liu model to study its property. As an 
example, taking m=6, the Liu PDF writes as: 



1 V5 1 V5(5 



X 



2^ 



' ttS + x^ tt (5 + x2)2 ^ ^ 

where the first term is a Cauchy law 6*1(0, v^, 0). As shown from the FiglHl the second term 
possesses a significant positive values at the middle and the medium parts of the probability 
law. The magnitude of these positive values are responsible of the height of Liu PDF and the 
departure from the stable law at the medium range values. The negative values of the second 
term of ( ITSj) explain the similarity of the asymptotic behavior of both laws observed in FigJTl 



For higher values of m, the same expansion holds with a Cauchy law S*i(0, \/m — 1, 0) as a 
first term but the second term becomes a sum of rational functions whose integral vanishes 
to ensure the normalization of the PDF. 



17 



10" 




FIG. 7: Semi-log representation of Liu PDF for m=4.5 and the stable law fit with a = 1.221 and 
7 = 1.3384. The error e = 0.066277. Solid line: Liu PDF model. Dashed line: Stable law. 

According to the probability theory, the modified Liu PDF belongs to the family of the 
Generalized Student or T-distribution whose expression writes as: 

'^^^''^ ^ wrm ((^^T^T^) ^^^^ 

where m is now called freedom degrees and a is a scale parameter. The characteristic function 
of the T-distribution is given by: 

latlf Ki^{\8X\) 
mt)= 2 ' '^^^^ll' (17) 

where Ki is the modified Bessel function. From ([2]) and (fTTI) . one can evaluate the difference 
between the stable and the T-distribution laws. As said in previous section, the stable laws 
decay as with a < 2. Here, it may be shown that the T-distribution laws have 
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FIG. 8: The corresponding figure of expression (jlSp . Tlie Caucliy law (first term) is on daslied 
line while solid line represents the second term with negative part. The zero value is drawn with 
dotted line 

the same asymptotic behavior as Ixl"™"^ for all m>0. Note that as a consequence, the T- 
distribution laws are more accurate than the stable laws to approximate the near-Gaussian 
situation. Another significant difference concerns the variance of the laws: for all m>2, the 
T-distribution laws have finite variance eiven by o"^ = while the stable laws do not 
possess finite variance except only for a = 2. From the practical viewpoint, this lack of 
variance may be of crucial importance when one has to compute statistical values as the 
mean square slope that is necessary to evaluate the radar backscatter cross section. 

Finally, from the expression ( JT6l) . it appears that the T-distribution does not include a 
skewness property while the stable laws are naturally able to express the skewness effect. 
Recalling the importance of the skewness effect on the dynamic of the surface waves, this 
point may have significant consequences on the accuracy of the slope distribution model. 
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V. CONCLUSION 



The main purpose of this study was to highhght the possible use of a-stable laws for water 
surface slopes representation. As said in the introduction of this paper, the fundamental 
argument on the use of a-stable laws stems first from experimental observations of the main 
property of the wind driven water surface waves: a possible fractal representation, related to 
a spectral power law, nonlinear characters of the wind wavefield resulting in an asymmetry 
and a modulational instability of the surface, and a randomness of the slope, such as these 
abrupt changes shown here in section 3. 

Considering these strongly disordered behaviors that may be found on slope variations, a- 
stable laws seem a priori to be potential candidates for this objective. Starting from this idea, 
the study propose to apply the stable laws to represent the probability distribution of the 
slopes of wind driven water surface from laboratory experiments. It is found that the stable 
laws are well adapted to model the random characteristic of the slopes data in different cases 
of experiment configurations. From the proposed cases, the four stable parameters behave 
in good agreement with the effects of different physical processes observed and expected in 
the experiments. In particular, the stability index traduces fairly the non-Gaussianity of the 
slopes in the same sense as the Stokes and/or the modulational instability effects. Also the 
skewness parameter reproduces well the asymmetry tendency of the slopes under various 
situations. However, at extreme experiment condition (large fetch and high wind speed), 
the stable laws fail to reproduce the behavior of the upper part of the slopes. 

From more theoretical consideration, comparisons with the Gram-Charlier and the Liu 
et al[l^ models are carried out. It appears that the stable law differs completely from the 
Gram-Charlier model as well as in its mathematical structure than from the characteristics 
of its results. Indeed, one can recall that the Gram-Charlier model at its usual truncation 
order (third or fourth order) is part of the near-Gaussian approximation models, while the 
stable law is totally non-Gaussian except only when the stability index reaches its maximum 
value. 

The comparison with the Liu model shows more similarities. Recalling that Liu model 
belongs to the T-distribution family, it is well established that both models are part of the 
so-called "heavy tail probability laws". They are both infinitely divisible probability laws 
but the T-distribution is not stable. This means that the sum of independent identically 
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T-distributed variables does not follow the same law unlike the stable law that it is always 
the case. 

From mathematical viewpoint such stability property lies onto fundamental characteristic 
based on the generalization of the central limit theorem (see Samorodnitsky and Taqqu 23l |) 
and traduces the fact that the a-stable laws constitute limit laws (or attractor laws) for all 
additive random processes. As a special case, we recall that for a = 2, from the central limit 
theorem, the Gaussian law is a limit law for all additive finite variance random processes. 
This is the profound justification of the popularity of the Gaussian model in many physical 
domains. In the realm of water surface waves, the random wavefield is usually considered as 
the sum of infinite number of random components. Owing to the weakness of the interactions 
between components, the corresponding motion of each component may be regarded as 
independent (see for example Phillips 1223). Consequently, under the generalisation of the 
central limit theorem, the choice of stable law is then theoretically justified with a clear 
mathematical basis. 

With regard to the practical viewpoint, the differences between the two laws are more 
important and have to be stressed here with respect to the study objective. These differences 
lie on four points concerning the closed form of the PDF, the finite variance value, the 
skewness representation and the interpretation of the physical processes through the PDF 
parameters. 

As said in previous section, the T-distribution has a closed form of the PDF for all values 
of its parameters . There are only three cases of stable distributions for which a closed form 
expression of the PDF is known. This lack of closed form is apparently an inconvenient 
property of the stable law. However, as shown in this study, there exists robust techniques 
to implement the integral form of the stable PDF (see Zolotarev 31 1, WeronjsQ], Nolan 19|). 
Other expansion techniques based on the Bergstrom series are also available to evaluate the 
stable PDF. 

For physical applications, an appeal to statistical quantities associated with the PDF 
model would be useful. For example, for the sea surface remote sensing problem, it is known 
that the mean and the variance of the slopes are of importance. This is historically related 
to the assumption of near-Gaussian state in the development. Recalling that T-distribution 
has finite variance for all freedom degrees greater than 2. T-distribution appears as a well 
adapted choice for this purpose. However, the T-distribution and the stable distribution 
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constitute fully non-Gaussian alternative and then fail to belong to the category of quasi- 
Gaussian solutions. In our opinion, an extension of the remote sensing or the hydrodynamical 
relationships in the sense as developed by Guerin 3] is required for a fully non-Gaussian 
solution. 

As far as non-Gaussian random water surface is concerned, the skewness is undoubtedly 
one of its most important properties. The stable distribution has asymmetric form that 
allows one to include naturally the skewness. This is not the case of the T-distribution for 
which no extension to asymmetric case is previewed. Empirical asymmetric formulation has 
to be found (see Liu et alfisl]). 

Before ending the paper, we focus on the interpretation of the physical processes through 
the probability distribution model. As we have shown in this study, the four stable param- 
eters are able to traduce correctly the effects of physical processes as the asymmetry and 
the kurtosis of the slopes or also its level of non-Gaussianity. Use of a-stable distribution 
instead of the traditional Gaussian or quasi- Gaussian distributions, combined to hydrody- 
namical laws will be of great help on understanding the physical behavior of the water wave 
slopes. In this sense, extension of the present work to the space representation of water wave 
slopes by the use of multivariate a-stable distribution has to be carried out and constitutes 
our prospective task. An other way to be explored is also the use of truncated stable law 
that allows one to obtain finite values of statistical moments. 
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